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Preface 


The factorial function known today under the name of Gamma function and the 
notation I(x), is due to the genius of mathematician and physicist Leonard 
Euler(1707-1783). It has been a subject of research since the XVIII century by some of 
the greatest mathematicians. 


To the following question: 
What is the function f of real variable x that has the propriety: f(n) =n! ? 


Euler found two solutions 


— the first solution consist in that f can be given as an infinite product of terms 


Pox) = He a. 

— and the second solution is that f can be represented as a definite integral 
b 

P(x), = J, g(x)dx. 


And since that, as we will see, and from the point of view of rigor, there has been no 
agreement on what should be the right definition. The reader will learn thorough 
these notes that there have been many points of view to define the celebrated 
function, and in a lot of references he will read its two definitions among statements 
like: the gamma function is the only function that verify such or such.... 
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Anyway, in modern terminology, Gamma function belongs to a large set of 
complicated functions used in many branches of modern science, formerly named 
transcendentals. Nowadays these functions are called special functions , very 
famous and useful functions like : theta function 0(z), Bessel functions , 
hypergeometric functions, Airy function, Hermite function... 

Clearly, by its origin, the Gamma function belongs to the world of mathematical 
analysis, but it is very important in many other mathematical contexts, either in pure 
or in applied mathematics. In pure mathematics for example, it appears in the theory 
of numbers, where it plays a preponderant role in the comprehension of the zeta 
functions linked to the distribution primes, and also in probability, and in theoretical 
physics... Regarding applied mathematics, it also appears in many domains like: 
Statistics, applied functional analysis, differential equations... 


The gamma function was first got in the context of real variables, then as many other 
important functions, has been considered for complex variables. In the case of 
gamma, the transition from IR to complex analysis in C was first accomplished by 
Gauss then used much by Riemann in his big work on the zeta function. So the image 
of the process around the transition from real to complex functions can be in the 


chart bellow: 
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The following notes, were originally annexed to my essay on the Riemann zeta 
function, and | have decided to expand and put it in an independent part, and | 
preserved the original title 


Mister « Lady Gamma ». 


| wanted these notes to contain some comprehensive information, both technical and 
historical, on this function presented in an easy manner to beginners. 


In fact the three main preoccupations of scientists that led to the calculus and other 
discoveries in analysis like our Gamma were: 


*Quadrature of the circle. 
*Finding areas and volumes. 
*Triangular numbers. 


So | will begin form the early works on the operation of interpolation which has led to 
the discovery of the gamma function, and the work done after the great Euler wrote 
its expression down. 


At the end of each part of these notes, | added some exercises. The exercises are 
related to the subject of each chapter; they are on calculus and its standard problems 
that turn around: convergence, calculation of definite and indefinite integrals, 
complex variables... These drills are meant to hint to the properties of this function in 
a way or in another, so | hope also that they help the reader to understand more this 
function. 
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Finite differences 





Before the definitive invention of infinitesimal calculus and its powerful methods, 
scientists used tricks to find values of functions by the operation of interpolation in 
the form of finite differences, a technic very fruitful and linked to Calculus as Georges 


Boole said in the beginning of his book 


In the following exposition of the calculus of finite differences, particular attention has been 
paid to the connection of its methods with those of the differential calculus, a connection 
which in some instances involves far more than a merely formal analogy. 


And further he defined it 


The calculus of finite differences may be strictly defined as the science which is occupied 
about the ratios of the simultaneous increments of quantities mutually dependent. The 
differential calculus is occupied about the limits to which such ratios approach as the 
increments are indefinitely diminished. 


What we known about this in common references is that interpolation existed in the 
work of Jacob Bernoulli (1655-1705), seen as the early typical example of 
interpolation of sequences but this mathematical philosophy of thinking about 
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numbers was not new at the time, there have been some earlier extensive work that 
treated interpolation . In fact, before Bernoulli was working on number combinations 
for probability theory, the language of interpolation was much used by the scientists 
of the XVI century. 


So the beginning of our story about the gamma function can be found in an activity 
much done by mathematicians, especially in the XVII century, an consisting in using 
clever formulas to generalize sequences of numbers, and this in general can be 
accomplished by two methods: 


1) Interpolation. The operation of finding values governed by a law within an interval 
between two real numbers a and b. 


2) Extrapolation. The operation of finding values, governed by a law, outside the 
interval bounded by two real numbers a and b. 


Simple examples of interpolation 


Some interpolation operations can very simple, for example to generate all odd 
numbers we may calculate the differences of the squares like this: 


During the early period known before steps to the invention of calculus, with for 
example the triangular numbers 


1, 3, 6, 10, 15, 21, 28, 36, 45... 


Which known to exist from Pascal’ triangle 


* This triangle have been discovered first by Muhammed Alkaraji(ca. 1050) then by Tshu Shi (1303)and last by 
Stifel (1544). 
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1 


11 

12 1 
13 3 1 
146 4 1 
15 101051 


16 15 201561 


The triangle shows three numbers: triangular numbers (in red), binomial number 
(horizontal numbers beginning by the 1s) and pyramidal numbers (in blew). The 
Pascal's law for difference 


(= CE) G 


Hariot and the others 


Many mathematicians of the XVII century like Napier (1-1) , Mercator(1-1) and Burgi( 
1552-1633)...were implied in practical calculations for navigation and other purposes. 
In fact the work done by those men can be seen today as a part of numerical analysis. 


Such mathematicians used interpolation in constructing: 
-Tables of sins for angles. 

-Tables of tangents. 

-Tables of logarithms, newly invented by Napier and Burgi 
-Tables of antilogarithms. 


English scientist Thomas Harriot’ (1560-1621) was an early English mathematician of 
the XVI century who influenced a lot of English mathematicians like Wallis and 
Newton.. He founded his work upon triangular numbers and to the laws governing 


* He has used the equal sign "==" similar to Record's " = " and we owe him the two signs <and >. 
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the binomials before Newton 


n(n — 1) n(n — 1)(n— 2) 
1.2 7 1.2.3 _ 


lon 
els 


The aim of Hariot was practical; he was implied in such mathematical activity mainly 
for the construction of navigational tables. 

His work consists in interpolating natural integers like those in arithmetic 
progressions by the study of constant finite differences 


Aty = AY = Vei — Ye 
and other higher differences 
A® 


as the example that precede. We can say here that as the derivate f’ gives us 
information about f, also the difference A” can led us to some valuable information 
on the sequence. 


Wallis 


Wallis (1616-1703) is also known for using interpolation and his work, we will see, his 
formula for : inspired Euler in his choosing the integral representation of the 


gamma function. Wallis was obsessed also by two problems: 
* of finding areas of figure (called quadrature). 


*and volumes(called cubature). 
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In his book Arithmetica infinitorum, he was concerned much of the search for the 
means of arithmetic progressions. He for example found that the means of the 


: . 1 
sequence of integers to the sum of last term is 5 


14+24+34+-4+n 1 


n+nt+n+t--+n 2 


an easy limit by the triangular number general formula. Some of the sequnces was 
hard, and he particularly asked help from the mathematicians of Oxford for the mean 
of the sequence n! 


1, 6, 30, 140, 630,... 


after futile attempts he concluded that such law governing the mean can neither be 
arithmetic nor geometric. In fact, he commented (p. 161) 


..and indeed | am inclined to believe (what from the beginning | suspected)that this ration we seek in 
such that it cannot be forced out in numbers according to any method of notation so far accepted, 
not even by surds... so that it seems necessary to introduce another method of explaining a ration of 


this kind, than by true numbers or even by the accepted means of surds. 


Wallis was interested in sums of sequences of powers 
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(1 — x8)! 


One problem especially occupied him and as well as generations of geometers 
before, it was the famous squaring the circle inherited from the ancient Greek. The 
search was the reason that led him to the quadrature of 


© = 


And thus to a ratio for the magic number 7t 





2X2X4X4X6X6X8X.. 1 
3x3x5x5x7x7x.. x4 


we see that this fraction has split n! in two parts : odd numbers in the denominator 
and even ones in the numerator . 


Trying to find such a relation in modern notation 
1 
1 
fa —x*)2dx =m 
—1 


He was implied in calculating definite integrals 
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and with the exchange t = cos(x) we can find 
1 
{a —_ ty\e-D/2¢-1/2 dy 
0 


and we will see that a similarin0 <x < 1to 


© 1 


{a —t)*t*"tdt = 


0 





sin(mx) 


a known identity for the gamma function. 


Newton 


The work of Wallis has influenced Newton(1643-1727) and led him to the binomial 
theorem. In fact before the posthumous work of Bernoulli, Newton himself inspired 
by Wallis and in 1665 , interpolated the very known binomial theorem for integers 


n 


(a+b)" = » (;,) ape 


k=1 


and he also went further to consider in place of the integer n the rational numbers x, 


but although (, 


operation of subtraction has a meaning 


; has by definition no meaning in this case when 7 is rational, the 
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x x(x —-—1 x(x -1)(x-2 
(1+ 1425422 pay FOE 


all this was published Newton's Methodus Differentialis in 1711. Euler later used the 
same formula when n tends to +00, to find the value of e 


1\" n n(n—-1)1 n(n—1)(n-2) 1 
e=(1+-) pg A ee 
n 12 n? 1.2.3 n? 


So we see here the appearance of the sequence n! in the terms’ denominators. 


Another problem seen also as an interpolation, is to find a polynomial P, (x) that 
passes by some finite number of points : suppose we have an unknown function or 
law f that takes the values 


Var Yn 


at the n points 
Mi uceg 


as polynomials are very good functions for their simplicity of calculation, we desire to 
find a polynomial that passes through all of them, and thus approximates the 
unknown f at those points. Two methods were discovered and have been used to 
solve this problem and both used finite differences, the first was by Newton(1676) : 


n 
(x — Xp) 
P(x) = Yo + > sry 
1 
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Where h = x — X. And the second by Lagrange(17) and his interpolation polynomial 
has the form 


(30 — X44) oe (KX — Xp) 


P(x) a Lu (x _ X;-1) — (x —_ Xn) 


There exists also in this domain interpolations that use rational and trigonometric 
functions in place of polynomials, developed later by Cauchy. The use of this technics 
is not at all old fashioned, it is still used in numerical analysis and since the XVII 
century many mathematicians used them in many important situations to solve 
problems , in 1860 Boole wrote a big book of more than 400 pages on this subject, 
still in use. 


Exercises 


1) make the differences of first ten consecutive cubes, what do you remark? 


n — ae — 


1) prove Pascal's law as a difference 
(= Gena) 


1) Show that 2” = Dip, } This number is the number of all possible combinations 


of n elements, and the number of all subsets of a finite set EF. 

1) Let (1,2), (—1,3) be two points in the plane R. Find the polynomial of degree 2 
that passes through these two points. 

1) prove the known triangular number formula: 
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n 
n 
f(n) =t, =) k=142434-4n= 541) 
1 


1) show that the limit as n — +00 


14+24+34+-4+n 1 


n+ntn+t-+n 2 


© 1) let the integral of the form 
B 


| sin? (x)cos?(x)dt 


a 





where at least one of the p,q _ is odd. Find suitable changes that lead to the integral 


B 
{a — t*)*tPdt 
a 
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The Bernoullts 





In the XVII century, mathematician Jacob Bernoulli (1655-1705) in jotting with 
numbers, was interested to interpolate sums of the form 


and find closed formula for them for given integer n and r , which means that we 
have only to put n = ny to get the searched sum . A special case of this formula is 
that of the known triangular numbers: 


n 
nN 
D KLSL42434-4n=S(t1) 


1 


where in this case r = 1, aS we have seen. 
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Another known formula when r = 2 


Ng(Np + 1)(2ny + 1 
124 224.32 boo ng? = Moo Te FY) 


if we think differently , Bernoulli was seeking alaw f(n) for the sequence 
1", (142°), (14+ 2" + 3°),... 


Thus doing this, he discovered the very famous Bernoulli numbers, and the law which 
Jacob Bernoulli looked for is : 


n 


d" =n, ; Ps 


1 


The ideas of Bernoulli on integers and combinations of numbers , with his formulas 
were published posthumously in a book entitled Ars Conjuctiri (1713), which means 
art of conjectures, the book dealt about probability. 


Goldbach 


During the following century, i.e. the XVIII century, the son of Jacob , Daniel 
Bernoulli (1700-1782) , and German Christian Goldbach (1700-1782) were occupied 
by the same problem and this operation was known to be called by the same name 
interpolation . So if we succeed in finding this such function we first :facilitate to give 
the values of this sequence for any natural numbers and avoid possible inextricable 
calculations, and second : we interpolate the sequence to other argumentss that are 
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not integers, and thus we may discover a new useful general function perhaps more 
important and promising than our sequence . 


Both the two mathematicians were searching this time for two functions: 


i) the first f interpolating the harmonic sequence u,, 





Which means: find a function sum 


(n) =1 +44 ee 
fn) = 2 n 


ii) the second g that give the hypergeometric sequence 


gn) = 1! 


For every natural integer n, so in this case the task was to interpolate the sequence( 
or progression) defined by the product 


n! = | | (n —k) 


O<k<n 


So in the case of the factorial , or more promising task, the formula that was originally 
valid only for integer numbers n , may be generalized and have values when the 
argument n ( or x in f(x) ) is not necessarily an integer : 


f(x) = x! 
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Or said otherwise, the task was to find a formula in function of n that easily gives the 
searched term of the sequence 


U, = n! 


without calculating it by other means, a task which Euler’s called in his prospective 
paper :” the most important problem in the subject of sequences””®. 





Exercises 
1) calculte 
10 10 10 
De DR De 
1 1 1 


1) in the spirit of Jacob Bernoulli, show that : 


hs 2 
yr _ ae + > 
2 
1 
1) it is said that Gauss when he was a pupil calculated the sum )°' k* =by writing : 
S, =14+24+3+4+::4+(n-1)+n 
S, =nt(n-1)++4+1 


complete the calculations to find the sum: S, = =(n +1). 


* The summary of the paper. 
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1) Show that 2” = Dip, (;,). This number is the number of all possible combinations 


, Le. the number of subset of a finite set. 


1) show that 
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When Euler entered the scene in 1728, he has just left Germany for Russia and heard 
of Daniel Bernoulli and Goldbach researches. 


He saw the problem very important, in fact in his coming paper(1831) we can read 
the following 





Euler, asking the same question: find a law f that gives the terms of the sequence 
Uu, =n! : 


Or find a function f such that : 


4 om ce 
He means "mathematicians". 
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He solved the problem in two papers, and to be we can retrace the events before the: 
*) Euler arrived at Russia in 1728 and began afterwards to think about the problem. 


*) He first found that the searched function f(n) = n! cannot be algebraic but 
transcendental. 


*) He corresponded with Goldbach about his thoughts on the 3™ of October 1729 
, and announced that he discovered a product for this function f(7). 


*) in another letter to Goldbach , dated on 8 of June 1830, he presented the gamma 
function as a definite integral. 


Algebraic and transcendental 


When Euler began tackling the problem(only 23 years old), he first observed : 





So he formulated the problem otherwise: he classified progressions in two categories 
: algebraic and transcendental. Euler supposed that sequence n! perhaps may be 
found by either algebraic or at least exponentially, then he changed his mind as we 
will see later on. 


i) algebraic : a progression that we can directly get the values from its general by 
finite number of operations like the triangular number t,, 
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n 
tpSAte2trstetn= Tmt) 


In this case for example, in order to find a term of the sequence )i k 
iF 3, 6, AO; cs 


it suffices to put value of n in =(n + 1). And there is no problem in saying that the 


function that gives the sum of the first n integers is 


foxy =" 


ii) Transcendental: a progression that is not algebraic , that has usually infinitely 
many terms. Which Euler named transcendental, he means a progression which 
transcends (just as everything in geometry which surpasses the powers of ordinary 
algebra is commonly transcendental (sic)) the easy nature of algebraic sequences and 
functions. 

He stated that the progression 


n! 


is difficult in nature and it is impossible to find an algebraic formula, i.e. finite like 
that giving any triangular number f,,, for it thus it is impossible by any algebraic 
expression. 


So , after the suggestion of Daniel Bernoulli, Euler communicated his early thoughts 
in a letter to Goldbach (1690-1764) dated on the 3™ of October, 1729, Euler 
announced that he has solved the problem by putting it in a product. 


He solved the other problem proposed by Bernoulli and Goldbach namely a closed 
formula for the sequence 
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and he found also a very interesting approximation 


n—+00 


1 1 
lim ( + 5 Ses ram log(n)) =V 0577 a. 





This is also named gamma and is very obscure number, we do not even if it is 
rational’. Curiously the two problems are related through the equality 


as we will see. 


The number n! as a product 


After this correspondence with Goldbach, Euler published an article (1730) in Latin 
untitled " De progressionibus transcendentibus seu quarum termini generales 
algebraice dari nequeunt " , which is in English® 


ON TRANSCENDENTAL PROGRESSIONS 
THAT IS, THOSE WHOSE GENERAL TERMS 
CANNOT BE GIVEN ALGEBRAICALLY 


and in another letter dated on 8 of January the same year 1730 he announced that he 
could present the function as an integral. 


° See the book by Julian Havil : Gamma : exploring Euler's constant, published in 2003. 
° Translated into English by Stacy G. Langton (1999). 
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Euler begins by saying that the progression u, = n! at infinity became like a 
geometric progression, this is because n! looks like 


n! = 2*1 x 3k2 x 5k3 x 74x |. 


He had the idea to introduce the geometric progression r” , he wrote the general 
term of u,, = n! as an infinite product: 


2 2 a ae AS 


u, = =— X ———_ X —— ——_ X —— 
1+n 2t+n 34+n 4A+n 


in fact this is simply equals by elimination of numbers to the product 


1 2 3 . 4. . 
= x x 7 
tin 1l+n 24+n 34+n 447n 














or 
k! 
L ==. 
kaa(k + n) 
So we have 
| . : x : x : x 1 
= i > XX —_————- ——_— — 
“0 1+0°2+0 3+0 4+0 
1! : : x : x : x 1 
= |! = ——_ K — XK — XK —— 
“4 +1 24+1 341 441 
2! : : x : x : 2 
= Z\ >= —— =X —_—_——_- —_——_— = 
“2 1+2°24+2 342 442 
3! : : x : x : 6 
= 53. SFX ——<—\_ ——_— —— = 
“3 Lie os Ge 4 
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The reader with knowledge of the infinite process, sees at once the problem of 
convergence of this infinite product. But Euler as it is known about him in other 
works, does not bother himself usually about the convergence of such infinite series 
and products, and may writes something like: 





(eee = 0 
2 3 7 


An infinite product in modern terminology and notation can be written 


+00 


| | n—- +00 
1 


In order that an infinite series 


os 


1 


converges, a necessary (but no sufficient) condition is that : 


lim u, = 0 
n—+0o 


For infinite product, also a necessary but not sufficient condition for convergence is 
that 


© Mohammed Fulano University of Constantine 
mohammedfulano@hotmail.com 


the product 


1 2 3 ¥ 4. . 
= x x 
tin 1l+n 2+n 34+n 447n 














is not defined when n is a negative integer. We, like Euler did, try to put some 
integers in the argument to find out the results: 

















n = 0,1,2,3. 
Then forn = 0Oandn=1 
47,22 _ 12h 
0! 140 £«°£:141 1 
forn=2 
ai 1.2%. 27+*.3- ; Z yoo ; 
= x Mae = 42a KO XK = 
14+2 2+2 13 24 3.5 
and alson = 3 
7 2° De 3s jae 3° 43 ’ 53 P 
= x ) i A Tel, Ce, Gree, Gee = 
14+3 2+3 4 45 34.6 44,7 


one must realize that the products are infinite, but Euler says that it always converge 
for real numbers(except for negative integers). Otherwise if n = 0, we have: 
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1 2 3 4 
—_ —-x —-x-—-—xX = 
12 3 4 

If n= 1 ,wehave: 
2 3 4 5 
XS KO KS KX i = 
2 3 4 5 


So at the time we consider the product as a finite product, and Euler added about 


this product 





So he consider as an interpolation and not only an expression for integers n € N but 


alsowhen x €Q. 





whose meaning will be apparent later on. 


Exercises 


1) are the following products converging: 
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University of Constantine 


i) ta +=) 

7 o i 

ii) T° — =) 

1) how can treat the infinite product 


L2° j2*.3% 32 ae abn 5 
= x ———_— x ————_- Xk ——__ X 
1+n 2+n 34+n 4+n 
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The appearance of the numbertt 


Euler went on by saying that the previous expression for n! i.e. the product 


1.2 2°4.3% 
x x 
1+x 2+x 








was not at all practical, or as he puts it "of any use for computing the terms..." : 
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but it is suitable to when n is a fraction so as Newton did with the binomial theorem. 


1 ; 
When Euler put the value x = a he was surprised to the appearance of the number 


7, i.e. he found that : 














Or as he wrote 





He then here considers it as an infinite product. He immediately recognized the 
formula, and he stated that he seen it before in the work English mathematician de 
Wallis’(1616-1703) by the hypergeometric series, the product found by Wallis in 
1665, namely 


2X2X4X4X6X6X8X.. 1 
3x3x5x5x7x7x.. x4 


or the same 


”in his book : Arithmetica Infinitorum. 
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+00 
2n2n TU 


1 @n-D@n-1) 4 


John Wallis as many other with circle and the number 7 before him French Francois 
Viéte(15-16) wrote in 1579 





discovered that the area of half-circle of formula y = V1 — x? onthe diameter is 1 
is equal to this product. 


ll 4 


Fig. The area under the curve: y = V1 — x? between —1 and 1 


Euler by saying that since S = mr“ then 
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So Wallis found 


S 2.4.4.6. 6.8.10 


(2r)2. 3.3.5.5.7.7.9.9 
When the diameter 2r = 1, we have 


— 2.2.4.4.6.6.8.10 


3.3.55.7.799 " ” 


and thus 


mt 2.4.4.6.6.8.10 
2  3.3.5.5.7.7.9.9 


to which Euler refers "corresponds with my series’. He knew then that the sought 
function that gives us n!, could not represented by an algebraic (by a finite number of 
operations ) nor by the function e~. 


In fact he 





To deal with the number y7z , it is probable [3] that Euler had the idea to use an 
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integral for the sequence n! , namely 


b 
vi= | foode 


So, in another letter to his friend Goldbach dated in 1730, he mentioned a possible 
relation with integration through : 


B 
n! = [ eed 


a 


And he called that relation with the quadrature of the circle. 


In fact he said: 





The idea of Euler by trying , search for a formula containing the factorial n!, he 
thought about some formulas which integral gives quadrature like 7 , as the formula 
used before by Wallis and other. Perhaps he can find 


B 
n! ~ [ nena — x9" 


Such that he chose 
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B 
| x"(1 — x)"dx 


a 


Where r can be an integer or not. He took this from the similarity with the known key 
integral 


B 
[x 2(1—x) 2dx=7 


a 








d d 
From the known la = arcsin(-) , We first make in = the change 


t=x- 7 and the result follows. 


He knew well the binomial theorem due to Newton (1665) that 


n 


(a+ b)" = » B gn pr 
n! 


k=1 


then he put 1 — x simply in place of a+ b 


(1—-—x)"= 





The developing differential dx 


n n(n — 1 n(n—-1)(n-2 
(1 —x)"dx = dx —-—xdx + Sr hae — Ss ae ++ 
1! 2! 3! 
multiplying by x” we find as him 
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n—-1 n(in—-1)(n-2 
a) ) ert ay ee Sa ) ray foe 


x™(1-—x)"dx = x"dx — A yrtldy 4 
1! 2! 3! 


Then by integrating 


: xrtt nx't4 n(n _ 1)x"+1 


r(1 _v\Nqy — = a See 
[xa iy ae = ee AGO). 40.G ao) 
0 


We know today that the integration of an infinite sum term by term is not always 
possible but Euler did not care. 


Euler did not chose 0 and 1 as the bounds of the definite integral arbitrarily, but he 
meant by a = 0 to eliminate any undesirable terms, and by 6 = 1 in counting. 


He then put in the result x = 1 and got 


1 7 n n(n — 1) —nn—i1j\(n—2) 7 


r+1 L@+2)'1.2.¢+3) 1.2.3.7 +4) 


so here when: 











1 0 O(n-1) 0(n-1)(n-2) 1 
n= 0 : —_— — ——————— = ee oe a 
r+1 = 1.(r+2) 1.2 .(r+3) 1.2.3 .(r+4) r+1 
ned: 2 4g 1) 
rt+-1. 0s 1.(r+2) 9 1.2.(r+3) = 11.2.3 .(r+4) (r+1)(r+2) 
2D 4 12 2.1 _ __2(2)(0) = 1.2.3 
7 " r41. 9 1.(r+2)) 1.2.(r+3) 1.2.3 .(r+4) ~ (r+1)(r+2)(r4+3) 
1 3 3.2 _ _3(2))(0) ee 1.2.3.4 





N=3 ica ig ig@s) 1asien )  ~ Gaemere 


He concluded that the definite integral 
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1 


[xa —x)"dx 


0 


generates a sequence of the general term 


n! 
(r+1)(r4+2)...(rtn+4+1) 








So 
1 
n! 
T(1 —x)"dx = 
fe Pee = Ga ane D 
0 
Or 
1 
n! 
r+n+1)]/x'(1-x)"dx = —_Mm 
( | ( ) (r+1)(r4+2)...(r4+n) 
0 
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The fraction £ and transcendental outputs. 
g 


So until now we have faced no problems when r =whole numbers, so to get 


transcendental numbers, what do we put in the place of r in 


n! 
(r+1)(r+2)...(r.tn+1) 





he replacedr_ by the fraction : to get 


n! _ftmtdg f f oan 
F+gft29)...f+ng) gt | #0 ion 


as (r+n+1) becomes 


f+(n+1)g 
g 


He says at this 





So the problem he faced was how to get rid of the term (r+ 1)(r+ 2)...(7. +n +1) 
and just find our formula 
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B 
n! = [ eed 


a 


which make appear the factorial n! naturally. So the solution is that in 


n! _ftmthg 


- f 
F+gf+2g9)...+ng) gt J #a-% - 





to put f =1andg = 0. 


But before that, Euler had to do other acrobatics and expand his analytic talent and 


JG 
Skills to get the result. So If we make the change x = y/+9 we find 


1 = 


=a g se 
dx yftg dy = ——ylt9dy 
f+g 


ft+g 
Then the formula involving the integral becomes 


tf DG 9 yh) Gh 
(f +g)(f + 2g)...(f tng) — grt | peal <P) xI+9 xft+9 dx 


but as 


and 


1 


ee ee AChE ee 
(ft+tgft2g)...ftng) ftg J a(t <P) dx 
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with also 


f+mtg_, 


lim Ff + Qf + 2g)... +ng) = lim rep 


g-0 g-0 





So the rest of the formula becomes 


n 


mee is, 
yim { (222?) ae 
J 


n! 





we find - . Thus Euler considered the derivation of this quantity by using the known 


Hospital rule to find 


d(1 — x‘) 


= —x'log(x) 
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and so 


1—x! 





= lim (—x‘logx) = —log(x) 


the appearance of the natural log , and if we return to the preceding formula, 











namely 
ee 
— 
we find 
1—x'\" 
: a a n 
um ( ; = (-log x) 
SO: 


IMatsxela=vanme 01 (=) am Wacie)) 


1 
a | Ctogxy dx 
0 





This formula got out of Euler's head is really amazing and is considered as one of the 
greatest discoveries in Analysis of the XVIII century. We can now try the three first 
integers into it to calculate the first three terms of the progression w,, = n!: 


1 1 
0! = | (-togx)' dx = | ax =| 
0 0 
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al 
1! = | -togx) dx =x —xlogx|5 =1 
0 


1 
2!= | logx.logxdx = 2 
0 


Amazingly, Euler generalize his formula for the function when the integers are 
n # 0,-1, —2,... 


n! 
r(x) = lim . 


n+0x(x +1)..(*« +n) " 


and also 


7 





x(x -1)!= 


sin(1x) 


called complements formula. 


Mathematicians of that epoch called the new gamma function a transcendental, 
which comes from the non-algebraic nature of the function I(x) ; we have seen a 
relation between the progression n! and transcendental functions. 


so the reflection formula: 


Exercises 
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1) find the area of by y= +v1— x? 
1) calculate: { log?(x)dx. (it was compute by Euler in his paper). 
1) by integration by parts, find 2! by Euler's formula 


if 
2S | -logxy dx 
0 





1) show that [ (=) = aa make an appropriate transform in J, (-logx)" dx to get 


+00 


| e~* dx 





1) calculate 
— 
1—-t? 
1) use it to calculate 
1 
[ia 
Vaz —t2 


1) Inthe, when a =O and 6 = 1, 


3 4 5 
(l—-x)"=1-—x4+=x%*-=x?4+- 


1! 2) 3! 

so find 

B 

4 

[a-xrax = 1-3+7t 

a 
1) calculate 
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a 
| a*—tdt 
0 


1) find @ and f such that we have: 





B 
_i _i 
[x 2(1—x) 2dx=T7 
a 
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At the same time an exactly in 1730 , Scottish mathematician James Stirling (1692- 
1770) was wondering about the same problem, but in a different manner. His book 
Methodus was occupied about numerical analysis. In his book he a chapter under the 


name "a treatise on summation and interpolation of infinite series”, he used finite 
differences to investigate, say about terms of a sequence (p.19): 


... When this have been found as Descartes has defined curves by algebraic equation: when 
these things have been obtained, problems about summation and interpolation and other 
matters 


He gave an approximation to the number n! known under the name Stirling formula 


Which is really good approximation of n! for big integers n. 


Using Stirling formula we can find that: 
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In(n!) ~ nln(n) 


and thus 


Hepergometric series 


In his Stirling was to a sum of an infinite series 


And he as sequences and wrote something like 


1 1-—n (1 —n)(2—-n) (1 —n)(2 —n)(3 —n) 
(oat), tere x(x +1)(%+2)(*+3) x(*4+1)(%4+2)(xX 4+ 3)(x + 4) i 


We see at once that Euler's is like a hypergeometric series. In fact if we take 
consecutive we find 
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Where the sequence has the form a” and when the is it is called a hypergeometric 


and for a function. 


Stirling 
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these series later saw another developments with Gauss and Kummer. 


Exercises 





© 1) give an approximation for 10000! 


1) show by using Stirling formula that: 
i) log(n 1) ~nlogn—-n 


ii)n" ~n! 
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A new expression 


Like the genius Euler, his gamma function gamma has occupied a lot of 
mathematicians after him. In his “Treatise of elliptic function and Eulerian integrals” 
(1826), Legendre introduced his reformulation of the new function, introduced the 
symbol I, and he named it "the integral of Euler’. 

With changing the variable ny = x inthe integral found by Euler, 


1 
n! = | Ctogxy dx 
0 


we find 
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Legendre wrote n! In his known notation ['(7) (third Greek letter) now known in all 
the literature, so we have: 


The Euler function is 


+00 


r(x) = i et 1dt 
0 


for x > 0. 





Of course the integral is seen as a limit 


It is not known why Legendre wrote x! instead of (x — 1)!, we can read in Edwards. 
p.8: 


... Unfortunately, Legendre subsequently introduced the notation I'(x) for I(x —1). 
Legendre's reasons for considering (n — 1)! Instead of n! are obscure (perhaps he felt it was 
more natural to have the first pole occur x = 0 rather than at x = —1) but, whatever the 
reason, this notation prevailed in France and, by the end of the nineteenth century, in the rest 
of the world as well’. 


For the convergence of I(x) we may write 


8 ° 
See Bronwein and Corless who have given some other quotes. 
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+00 1 +00 
| e't*1dt = [etetac+ | e*t*ldt 
0 1 


0 


And we have 





1 1 
1 
O< [ etertat < | t*!dt = ~ 
0 1 
Andalson >x-—1 
+0o +0o 
0 <| e*t*ldt <| e§t"dt < +00 


1 


nN 
The latter is calculated by integration by parts and knowing that lim;_,,.. — = 0. 


It is easy to show that ['(n) = nI'(n), in fact by a integrating by parts the integral 


T(x) = [ etetat 
2 
we find 
x 
T(x) = —e* + | e't*ldt 
2 
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which gives 


(n+ 1) = nI'(n) 


or more generally the known relation that gives n! = n(n — 1)... 1: 


So in real variables 





r(x+1) =2T (x) 





This formula is known under the name " the functional equation of [' "n dit is know 
to Euler and it aprs in the functional equation of the zeta equation. 


lt enable us to prolonger the function I for negative values , thus 


rix+1 

rey = Pet D 

X 
for example give the values 

r'(0,5) 

r(—0,5) = —— 
(—0,5) 05 
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For x < Othe function I’ looks like the functions 


sin(x + n) 


where € R, taken from the formula 


7 





r(z)V¥(4—-—z)= 


sin(1z) 


Called complement formula. Legendre found the formula de multiplication 





r (2) (224) - 2 re 


Known as Legendre's duplication formula, as we shall see a special case of Gauss’ 
multiplication formula 


The Beta function 


Legendre named I(x) : Eulerian integral of second specie. And the first specie is beta 
function 


1 
B(x, y) = fa — t)’ 1t1dt 
0 


Which converges for x > 0, y > 0, and it is easy to see that 


B(x, y) = BW, x) 
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The other Eulerian introduced by Legendre was beta function is linked to I’ through 
the formula 


Tory) 


B(x,y) nea) 





Binet 


Two formulas attributed to another French mathematician, namely Jacques 
Binet(1786-1856) , the first is 








+00 
logl'(z) = (z-~)I +~log2 +{ (6 i ) a 
ogI'(z) = (z >logz—z+ logan a ee x 
0 
and the second 
00 x 
1 1 arctan(—) 
logI'(z) = (z- 5 logz —Zz+ 5 log2n PZ | aa 
0 
Exercises 
1) using the integral find (7). Compare the result with 
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1 
log(x) 





1) using show that 


transcendantal. 
1\_ 7 (3\ _ vz 
1) show that (=) =[T (=) = 
1) show that B(x, y) = B(y,x). 
1) prove that the Beta function converges when: x > 0,y > 0. 


1) show that 





B(x, y) = BU, x) 


1) using B(x, y) integral definition , show that 


ro) 
BUX, = —— 
(x, y) T@+y) 
Hint: use first [(x)I(y). 
1) show that 
1 
B(x,1) =- 
X 
1) show that 
TU 
BGS). = —— 
sin(7x) 
1) show that 
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As a law, many important functions in analysis are really interesting and give heir 


hidden pearls when they are considered for complex variable z. we will learn here 
that it was Gauss who considered for the first time the function I’ for complex 
arguments. 


Complex numbers for [ 


In 1823, and in his work on infinite series, German Gauss introduced the new 
notation II(x) in place of (x + 1)! antithetical to Legendre 


By the fact that 


M(z-—1) =T(z) 


or for R(z) >-1 
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M(z) =f. e-‘t’dt 


this definition has been the one used by Riemann in his celebrated paper on the zeta 
function, one of the rare writers of mathematics who used it was Harold Edwards in 
his analysis of the Riemann paper(1859) p.8 , where he remarked: 


Gauss’ original notation appears to me to be much more natural and Riemann's use 
of it gives me a welcome opportunity to reintroduce it. 


Gauss preferred also to define the gamma function as a product 


n!' 
Ii(z) = lim 


n—+0 (z+ 1) ~ (z+n) wer ay 


Which may be defined for all complex values except the negative integers numbers 
(which make the denominator vanish) this formula we have seen is known to Euler. 
Gauss wrote to his friend German Bessel(1748-1846) : 


But if one does not want... countless fallacies and paradoxes and contradictions to be 
exposed, 1.2.3 ...x must not be used as the definition of (x), since such a definition has a 
precise meaning only when x is an integer; rather, one must start with a definition of greater 
generality , applicable even to imaginary values of x, of which that one occurs as a special 
case. | have chosen the following... when k becomes infinite’. 


and some mathematicians prefer Gauss definition , and it was Liouville’s []. 


Gauss also discovered the multiplication formula 


7 R. Remmert: Classical topics in complex function theory (p.34). 
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r(=)\r(—) te) oe (x) 


When p = 2 it is Legendre relation seen before. 


The analycity of ['(z) 


It has been proved that the function of complex variables: 





+00 


l'(z) =| e~'t2-ldt 


0 


is an analytic function except for simple poles at the points: z = —n. 


we the theorem 


Theorem 


The function I is defined and analytic over all complex numbers C , except for simple 


poles at z = —n with residues 





We can now state the known formulas of the gamma function in the complex case: 


i) '(z+1) = z2T(z) 


7 


ii) (z)rQ —z) = 





sin(mZ) 
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| 
ii) T(z) = lim ep Te (n + 1)” 
i (2) = 2r (2) Ea? 
ii) '(z)T (z + ~) aT (z + “—) = (20) V/2nG/2)-n27 (nz) 


And in Gauss’ notation 


i) T(z) = zll(z — 1) 


ZU 


ii) W(z)M(—z) = 





sin(mZ) 
" , nt 
= tet) @en rt) 





v) H(z) = n7T] (=) I] (=) Tl (=) ant which gives the Legendre formula for 


n (21)” 


n=2 


iv)II(z) = 271] (=) nn? 


Exercises 
1) prove that 
= 1-2 Z = 
r{" (n + 1) [Ja " 2) 104 : yz 
zt+n 7 n n 
1 1 
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1) show that T (=) = Vr =Tl (- =) 


1) calculate by the definition of Gauss that I1(0) = 1, what is the value of ['(0)? 
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Weserstrass 





one of the preoccupation of German mathematician Karl Weierstrass (1815-1898) 
was to write a function in function of is zeros or roots. The function I(x) has no zeros 
, but as it is not defined at the points = —2n, and by theorem 


Al 


If lim, f(z) = © then lim, a 7 = 


And the function = should have zeros at the point a. Likewise the function a 


Z=-—2n, 


As the function is analytic like a polynomial the idea was to seek a possible result that 
assure the representation s a like polynomials themselves, as a product in function of 
the roots zeros 


| |a — ZeTOS,) 
k=0 


which gives: 
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Theorem of Weierstrass(18): 


1 YZ ; (1 a “) —Z/k 
= Ze oa 0&4 
I'(z) | | k 





Remmert refers the above formula to German Oskar Schlomilch (1823-1901) in 1843 
and to a certain Newman(1848), whom | failed to recognize’’ . He added : " It has 
become customary to call ...the Weierstrass product but it does not appear in this 
form in his work...". 


However Weierstrass in one of his articles wrote 





1 — x i n+1 
— | [4D ects 
r(x) | | pe" 


where y is a number we have seen before, and it is called Euler-Mascheroni constant 
or simply the constant gamma , 


1 
y = lim ) = — logn) 
n— oo nN 
it has been calculated by Euler to decimals 
y = 0.577. 


We do not even know if itis rational but the relation with function I is also 


° It is surely neither German Franz E. Newman(1798-1895) nor his son Karl Newman(1832-1925). 
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(1) =-y 


Since the function ['(z) has poles z = —2n, then we infer that — must have the 


zeros Z = —2n, and in another formula due to Weierstrass. 





Curiously, Riemann wrote once 


A-Ag=+oo 
n= 11 | | pee Ge by 
we shee ra—n+a) T(1—a) 
A-Ay=—-0 


Where 0< a) <1and 


Holder 


in 1887 , German mathematician H6lder (1859-1937) proved : 


Theorem (1887): 


The function I’ cannot be a solution to any algebraic differential equation of the form 





© Mohammed Fulano University of Constantine 
mohammedfulano@hotmail.com 








He showed that any solution to the presupposed algebraic differential equation 


would not verify the typical functional equation of I’, namely 
T(z +1) = zI(z) 


Holder's result means that the function I is a transcendental function. This result is 
similar to the Hilbert result (18) on the function ¢, what shows that the two 
functions are of more complicated nature then algebraic function. We can define an 
algebraic function as that an equation like 


p+pt+ta=0 


And an transcendental function a function that does not an equation. Transcendental 
functions like cos(x) , sin(X),... 

some other special functions like Bessel functions , hypergeometric functions, Airy 
function, Hermite function.... Are solutions of linear differential equations. 


Kummer 


in 1887 , German mathematician Kummer (1859-188) proved that the series , and 
integral representation : 


F(a,b;c;1) = l(c — a)I'(c — d) 


this is a hypergeometric series by Stirling and studied by Euler and especially Gauss. 
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Hadamard 


In 1894 , French mathematician Jacques Hadamard found entire function f that 
coincide with I when x = n. 


at-\ek-leat-1 ae me) manlel tom @Rstoy-9): 


If x is real then 


2—xX 


—) 


Z 1 Z neues - 
f(x) = T(x) * dx og| (ZV ( 


such that f(n) = nl. 





Stieltjes 


In 1894 , Belgium mathematician Stieltjes (185-195) wrote an interesting article on 
the gamma function giving a useful estimation of logI'(z) 


Theorem 


+00 


1 
log! (z) = (z + 5 | logz—z+log2 — | rat 
0 


It is a very interesting approximation in practice, which we can roughly write 


logI'(z) ~ zlogz 
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this formula is much used in the calculation of the roots of the Riemann zeta 


function. 


The function 


logI'(z) 
whose derivative is named the bigamma function is: 


dlog!(z) _ I''(z) 
dz '(z) 


we curiously find: 


_ ra 
~ (1) 





Already seen relation relation between the two gammas. 


Exercises 


1) write the following in their zeros € C: 

i)x* +2x-1 

ii) x* + 2ix +i 

iii) x*° —2x+1 

iv) xt +2x7°+x4+2 

1) write ['’(1) in the integral form. 

1) in the product for [ take logI’ and expand it. 
1) show that I(x) is positive for all reals x > 0. 
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(Zz) nz, z-1 ue 
1) show that a 2°m* tang ( ) 
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Another definition founded on another point of view, is to consider the gamma 
function as a very special convex function. This "philosophy" was introduced by two 
Danish mathematicians: Herald Bohr’ (1887-1951) and Johannes Mollerup(1872- 
1937) in their book on Mathematical Analysis** , it was to define the gamma function 
by logarithmic convexity”. 


The idea of Bohr and Mollerup is that is logarithmically convex: any function that 
looks like gamma i.e. verifies the functional equation 


T(x) = xP(x) 


and is not logarithmically convex is a pseudo gamma function. 


have showed the following fact about I 


** Brother of the known Niels Bohr(1885-1962) one the fathers of Quantum Mechanics. 
* Tome Ill p.149-164. This book was written in Danish and published in 1922 (Copenhagen). 
** See Remmert(p. 44). 
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Bohr-Mollerup theorem (1922)8 


The function I is convex. 





And also 


Bohr-Mollerup theorem (1922) 


The function I is the unique solution to the equation [T(x + s) = xI(x). 





Convex functions 


Convexity is a propriety much used in Analysis and Geometry : in geometry for 
example , a convex set is a set in which any line segment between two of its points is 
still inside the set , geometrical figures that are convex : disk, triangle , squares, 
rectangles, ... a circle is not convex and any polygon is not necessarily convex. 


In analysis we can generalize this concept of convexity to functions by introducing the 


definition : 
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Or equivalently: the function f is convex if the set’’ of points in the plan above the 
f graph is convex. 


Fig. f(x) = x* is a convex function. 


if f:1 C R—R isa function two times continually differentiable over an interval 
then f is convex iff 


f(x) 2 0 


We also have these properties 

* The sum of convex functions is also convex. 

* The limit of a convergent sequence is a convex function 

* a convergent infinite series of convex terms has a convex sum. 


A function is concave 


“* This set is usually named the epigraph of f. 
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Fig. g(x) = —x? is not convex, but a concave function. 


The Gamma function verifies this propriety of convexity , for from Leibnitz's law we 
can write 


+00 


r(x) =| e~*(logt)*t*-*dt 


0 


and as obviously e~‘(logt)*t*—* > 0, the propriety follows. Besides, we can at once 
infer that from its graph in the real plane. 


Artin's Book 
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Artin began his book (1931) 


| have written this monograph with the hope of filling in a certain gap which has often been 
felt to exist in the mathematical literature. Despite the importance of the gamma function in 
many different parts of mathematics, calculus books often treat this function in a very 
sketchy and complicated fashion. | feel that this monograph will help to show that the 
gamma function can be thought as an elementary function, and that all of its basis 
properties can be established using elementary methods of the calculus 


Artin chose to treat the function I(x) by using the integral and the preceding theory 
of convex functions, he added: 


../| have chosen the integral as my original definition of the gamma function because this 
approach saves us the trouble of proving the convergence of Gauss’ product. Any other 
analytic expression having the characteristic proprieties of the gamma function could just as 
well have been used. The whole theory will then be deduced using the concept of log 
convexity. This method comes from Bohr and Mollerup. 


He a theorem from the theorem of Rolle. 


3 


the function f is said to be convex if, and only if f has monotony increasing one sided 


derivation. 





Then comes the previous result in the case of numerical function as a consequence: 
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oli ay : 


Let fa twice differentiable function, then f is convex if, and only if f(x) = 0 for all 


the values x of the interval. 





He then defined any weakly convex function as a function verifying the inequality 








1 
F() <5¢@+ro» 


and proved : 
: 


Let f be a is convex iff it is continuous and weakly convex. 





So we have: 
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Lemma 


If the function f is logarithmically convex then f (x) is convex. 





The inverse is false. 


Artin a un important theorem know under the name Artin-Bohr- Mollerup 





The function gamma is the only convex function that verify : 





Bourbaki1 


This point of view has been followed also by the Bourbaki school to characterize the 
gamma function . It is a generalist philosophy proper to the French group of 
mathematicians. 


The point they (and as did Artin p.13) is that defining gamma by the functional 
equation 


(x) = xT(x) 


We can infer infinitely any function that verify it such as 


rixt+n)=(*+n—-1)(x+n-2)..(*« +1) (x) 
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and 


I(x) (x) 


et ae ee eee 


In their book , volume IV (1976) , after they defined convex functions (p. 22) they 
defined the gamma function as 


There exists one and only one convex function f defined on |0, +00| that satisfies the 
equation 


f(x + 1) — f(x) = logx 


And take the value O for x = 1. 





Exercises 


1) show that 


. n! 7 
Me eee za +1). en) * 
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converges for all z except: 0,—1,—2,-—3... 


1) show that [(1) = 1 


1) show that the function (x) has no zeros in the interval ]0, +00[ by using 


+00 


fi ertemia 


0 


1) show that 
+00 
| e't™dt < +00 
1 


1) Show that lim,,4. — = (). 


1) show if the following functions are convex 


a) f(x) = —x* 
b) f(x) = x° 


1) show if f(x) is convex then e/™ is log convex. 
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